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ABSTRACT 

The present status of the use of two-particle intensity interferometry as a diagnostic tool to study 
the space-time dynamics of intermediate energy heavy ion collisions is examined. Calculations for 
the two-proton and two-pion correlation functions are presented and compared to experiment. The 
calculations are based on the nuclear Boltzmann-Uehling-Uhlenbeck transport theory. 



KEYWORDS 

Heavy ion transport theory; two-particle interferometry; two-pion correlation function; two-proton 
correlation function; pion shadowing; nuclear equation of state; in-medium nucleon-nucleon cross 
sections. 



INTRODUCTION 

Probably the premier goal in performing heavy ion collisions at intermediate and high beam energies 
is to obtain information on the nuclear matter phase diagram, on how nuclear matter behaves under 
compression and heating, and on possible phase transitions in nuclear matter. Trivially, the time 
development of the process with which we hope to compress and heat nuclei, the collision of two 
heavy ions, is not directly observable to us. To obtain this information we are exclusively dependent 
on the particles emitted during the course of the heavy ion reaction. Some of these particles are 
created in the reaction process (e.g. photons, pions, kaons, etas, anti-proton, . . .). Other particles 
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have already been present in the initial stages of the reaction (protons and neutrons) and only have 
their momenta and energies modified during the collision. 

All of these particles are observed a long time after the reaction has taken place. We are only able 
to observe the final four-momenta of these particles, but not the time development that led to these 
final states. Furthermore, we are not able to directly measure the location from which these particles 
were emitted. Thus we are confronted with a scenario in which we can only observe information 
which is integrated over space and time. It is our task to unfold this information to reconstruct the 
history of the heavy ion reaction and to possibly extract the physically relevant parameters of the 
nuclear equation of state. 

In this sense heavy ion reaction physics is like forensic medicine, in which minute tell-tale signs and 
pieces of evidence left behind are collected in order to reconstruct a crime. 

The last decade has seen an extensive investigation of inclusive single particle spectra from heavy 
ion collisions. The basic idea (Stock et al, 1982) is that compression of nuclear matter requires 
energy, which is then not available for the production of secondary particles. By first comparing 
to intranuclear cascade simulations (Stock et al, 1982) and then to nuclear transport calculations 
including mean field effects (Bertsch et al., 1984; Aichelin and Ko, 1985) it first seemed that there 
was a clear signal for the value of the compressibility of nuclear matter. However, further studies 
including momentum dependent nuclear mean field potentials (Aichelin et al, 1985; Gale et al, 1987; 
Gale, 1987), medium effects on the elementary scattering process (Cugnon et al., 1987), and improved 
elementary hadron-hadron cross sections (Li and Bauer, 1991, 1991a; Danielewicz and Bertsch, 1991) 
have shown that the connection between nuclear compressibility and produced particle spectra is not 
as striking as previously thought. The present status of the theory of particle production (Cassing 
et al, 1990) indicates that there is at most a factor of 2 difference in yield of produced particles if 
one changes the nuclear compressibility constant, k, between values of 200 and 380 MeV. This small 
sensitivity to the nuclear compressibility is due to the fact that the amount of compressional energy 
stored during the heavy ion collisions is roughly the same for all values of k. Instead of storing more 
energy in compression, nuclei rather adjust the maximum density reached in heavy ion collisions. 

The basic problem with examining single particle inclusive cross sections is again due to the fact that 
we are only observing a time and space integrated quantity. This problem is somewhat less severe for 
high energy photons, because they are produced very early in the reaction (Bauer et al, 1986) and 
have very little final state interaction. But pion and to some extend even kaon spectra are strongly 
modified by final state interactions. In addition, at beam energies of around 1 GeV A, the production 
process of these particles can in some cases be dominated by collisions of nuclear resonances in the 
nuclear medium, i.e. by processes for which we do not know the elementary hadron-hadron cross 
sections from experimental data. Thus inclusive particle production cross sections carry only limited 
information on the space-time development of heavy ion reactions. 

One can then ask if two-particle correlations carry more information on the space-time history of 
the heavy ion reaction process. The motivation for this is that by observing single particles in the 
final state one integrates the dynamical information over their entire world line, whereas by using 
two-particle correlations one has two independent world lines, the history of which one integrates 
over. Since the relative wave function of the two particles is modified due to their interaction at the 
point at which the world lines cross or come close enough (if they do), two-particle correlations can 
serve as a possible source of information on the space time development of heavy ion reactions. 

In this paper I will focus on the utilization of two-particle correlations at small relative momenta, at 
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which intensity interferometry is possible. 



Intensity interferometry was introduced by Hanbury Brown and Twiss (1954, 1956, 1956a) as a 
technique for astronomical distance measurement. They recorded the two-photon correlation function 
for incoming coincident photons as a function of their relative momentum. This correlation function 
can be written as: 



R(h,k 2 ) 



(ni)(n 2 > 



where (n^) is the probability of detecting two coincident photons of wavenumber k% and k 2 in 
detectors 1 and 2, and (rii) is the probability of detecting a photon of momentum k% in detector % 
(i = 1,2). Equation [I] contains only count rates, which are proportional to the absolute squares of 
the amplitudes. As a consequence, HBT interferometry is insensitive to phase shifts introduced by 
atmospheric disturbances. It can be used with very large base lines and delivers superior resolution. 
This was first shown (Hanbury Brown and Twiss, 1956a) by measuring the angular diameter of Sirius. 

The physical basis of the HBT effect is that two photons have a non-zero correlation function due 
to the symmetrization of their wave functions, a consequence of the quantum statistics for identical 
particles. 

A similar technique can also be used for source size determinations in subatomic physics. This was 
first utilized by Goldhaber et al. (1960) in studying angular distributions of pions in pp annihilation 
processes. They found that the emission probability of coincident identical pions is strongly affected 
by their Bose-Einstein statistics, which causes an enhancement of the correlation function at zero 
relative momentum, q — 0. The width of the maximum at q = depends on the radius of the 
interaction volume (Goldhaber et al, 1960) and also on the life-time of the emitting source (Shuryak, 
1973). 

In recent years, two-pion intensity interferometry has been strongly pursued at ultra-relativistic 
energies, where the interest is on the interplay between source dynamics and final state interaction 
(Pratt, 1984; Kohlemainen and Gyulassy, 1986) and pion correlations from an exploding source 
(Pratt, 1986). Intensity interferometry derives its main attraction, however, from the prospect of its 
possible use as a diagnostic tool for the formation of quark-gluon plasma (Pratt, 1986; Bertsch et 
al, 1988). 

Intensity interferometry is not restricted to bosons, but can also be applied to Fermions. Koonin 
(1977) proposed to use two-proton intensity interferometry to obtain 'pictures' of heavy ion collisions. 
The advantage of using protons as a probe lies in the fact that they are already present in the colliding 
nuclei and can be liberated relatively easily. In contrast, to create a pair of pions one has to spend an 
energy 

-^min — ^ TTl-ff ~ 280 MeV in the center of mass of the generating system. Therefore, protons can 
be used as a probe at much lower energies. In addition, the two-proton relative wave function contains 
the prominent 2 He-'resonance', which leads to enhanced sensitivity of the correlation function to the 
source size. Finally, protons are easy to detect with the required resolution. 

Recent progress has been centered around the theoretical computation of two-proton correlation 
functions from nuclear transport theory (Gong et al, 1990, 1991, 1991a; Zhu et al., 1991, Bauer et a/., 
1992, 1992a, 1992b). In this framework, it is possible to understand the dependence of the correlation 
functions on the spatial dimension of the emitting source, on the momentum distribution of particles 
in the source, and on the time development of the system emitting the particles. Comparisons of this 
theory to experimental data have established two-proton intensity interferometry as a quantitative 
tool to study heavy ion reaction dynamics. 
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A summary of the present status of the field as well as further references can be found in Boal 
et al. (1990) and Bauer et al. (1992). Here I primarily focus on the use of two proton and two 
pion intensity interferometry as a diagnostic tool for reconstruction of the space-time history of 
intermediate energy heavy ion reactions and with it the determination of limits on the parameters 
of interactions of nucleons in dense nuclear matter. 



HEAVY ION TRANSPORT THEORY FOR INTERMEDIATE ENERGIES 

During the last few years, several groups have developed a nuclear transport theory for intermediate 
energy (20 MeV < E^, C£Lm /A < 2 GeV) heavy ion reactions (Bertsch et al., 1984; Kruse et al., 
1985, 1985a; Gregoire et al., 1985, 1987; Bauer et al., 1986; Li and Bauer, 1991; Danielewicz and 
Bertsch, 1991). This transport theory describes the time evolution of the nuclear one-body Wigner 
distribution f(f,p,t) under the influence of the nuclear mean field and individual nucleon-nucleon 
collisions via the Boltzmann-Uehling-Uhlenbeck (BUU) equation 

r\ — » 

-/(f,p,t) + P-W r f(f,p,t)-W r UW p f(f,p,t) (2) 
at m 
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6 (2m ^ + g 2 ~ 9i' - <?2')J • + Q2- qv - ql') ■ ^ 
7(f, q v ,t) f(f, q 2 ,,t) (l - f(f, p, t)) (l - f(f, q 2 , t)) 
-f(f, p, t) f(f, q 2 , t) (l - f(f, q v ,t)) (l - f(f, q 2 ,,t))y 
Here U is the mean field potential, which is commonly parametrized as a density functional 

where the constants a, (3, and a are determined from the choice of a nuclear compressibility, k, and 
from the conditions E/A = —15.75 MeV and d(E/A)/dp = at p = p . Popular choices are (Bertsch 
et al, 1984) a = 7/6 => a = -358.7 MeV, (3 = 304.6 MeV ('soft', « = 200 MeV), a = 2 =>• a = 
-124.4 MeV, (3 = 70.3 MeV ('stiff', re = 380 MeV), or (Bauer et al, 1986) a = 4/3 =>• a = -218.1 
MeV, (3 = 164.0 MeV ('medium', re = 235 MeV). da/dQ is the energy dependent in-medium nucleon- 
nucleon scattering cross section, which is taken from free space elementary cross section data where 
available, and from detailed balance and isospin symmetry arguments otherwise. 

All present solution schemes for the above equation employ the test particle method (Wong, 1982), 
in which the entire phase space is divided into small cells, whose equations of motion are first order 
differential equations in time, 

i = l,...,(A t + A p )H, (6) 

where A t and A p are the target and projectile masses, respectively, and M is the number of test 
particles per nucleon (usually taken > 100 to reduce artificially generated numerical fluctuations). 
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Fig. 1: Time evolution of the coordinate space density p(x,0,z) in the 
reaction plane (in units of po) for the reaction 93 Nb + 93 Nb at a 
beam energy per nucleon of 60 MeV and impact parameter as 
calculated from the solution of equations [| and ||. 
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The term C(pi) represents the solution of the collision integral via an intranuclear cascade (Cugnon 
et ai, 1981, 1982; Cugnon and Lemaire, 1988) for the test particles. The test particle collisions 
respect the Pauli exclusion principle due to the presence of the factors (1 — /), which are numerically 
implemented via a Monte Carlo rejection method. In our particular numerical realization the values 
of f(r,p, £) are stored in a six-dimensional lattice so that the computation of the factors (l—f(r,p,t)) 
only requires the call of 2 6 lattice elements for a six-dimensional interpolation (Bauer, 1988). 

Further details of the numerics and physics of heavy ion transport theory can be obtained from the 
review articles of Stocker and Greiner (1986), Bertsch and Das Gupta (1988), Schuck et al. (1989), 
Cassing and Mosel (1990), and Wang et al. (1991). 

The solution of equations £| and || provide us with a time evolution of the single particle phase space 
density f{r,p, t) and thus yield (within the model and on a semi-classical level) complete information 
on the heavy ion reaction history. At present this is the best and probably only valid way to visualize 
the nucleus-nucleus collision process. By comparing the predictions of the model to experimental 
data one can iteratively check the approximations entering the simulations and refine them. 

For illustration we display in fig. 1 the time evolution of the reaction 93 Nb + 93 Nb at a beam energy 
per nucleon of 60 MeV and impact parameter. Displayed is the nucleon density (in units of po 
in the reaction plane (y=0). For this simulation (Bauer et al, 1992c) a number of A/"=1000 test 
particles per nucleon was used. The beam direction is along the z-axis. We can see the two nuclei 
touching (£=10 fm/c) and interpenetrating and compressing each other (£=40 fm/c), an expansion 
phase (£=70-100 fm/c), and finally the formation of a ring-shaped remnant. 

It should be noted in this context, however, that the model does not contain the formation of 
complex intermediate mass fragments other than those resulting from the decay of projectile and 
target remnants. A large number of groups presently work on the inclusion of fluctuations into the 
dynamics in order to overcome this difficulty (Bauer et ai, 1987; Ayik and Gregoire, 1988, 1990; 
Randrup and Remaud, 1990; Burgio et al, 1991, 1992; Bonasera et al, 1992, 1992a; Randrup, 1992; 
Reinhard and Suraud, 1992; Ayik et ai, 1992; Benhassine et ai, 1992). The verdict on the success 
of these attempts to describe nuclear fragmentation is still out, and it may well be that a true N- 
body quantum theory is needed. Present attempts to approximate this theory include transport 
models with dynamical production of A < 3 fragments (Danielewicz and Bertsch, 1991; Danielewicz 
and Pan, 1992) and on so-called 'Quantum Molecular Dynamics', semiclassical N-body simulations 
(Aichelin and Stocker, 1986; Aichelin, 1991; Ono et al, 1992; Peilert et al, 1992). 



SINGLE PARTICLE PROTON SPECTRA 

Since the solution of Equation ^| represents the time evolution of the single particle distribution 
function f(f,p,t), it is possible (in the limit of £ — > oo) to predict all single particle observables such 
as proton spectra (Aichelin and Bertsch, 1986; Bauer, 1987a) in this theory. 

As one example we show in fig. 2 the comparison of theoretical calculations and experimental data 
for the single particle inclusive proton production cross sections. The calculations are represented 
by the histograms and are the results of the calculation of 40000 different events at random impact 
parameters. The data (Fox et ai, 1986; Chitwood et ai, 1986) are represented by the circles. As one 
can see, the overall normalization as well as the energy and angle dependence of the experimental data 
are well reproduced. There is a maximum disagreement of about a factor of 2 between calculation 
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Fig. 2: Inclusive single proton energy spectra for the reactions 16 O + 12 C 
at E/A = 25 MeV and 12 O + 12 C at E/A = 40 MeV. The calcu- 
lations (Bauer, 1987a) are represented by histograms and repre- 
sent the result of the simulation of 40000 reactions at randomly 
distributed impact parameters. The data (Fox et ai, 1986; Chit- 
wood et ai, 1986) are represented by the circles. 

and measurement. Thus we can be confident that the phase space distribution function integrated 
over coordinate space, / d 3 r f(f,p,t), reproduces the experimental observables in the limit t — > oo. 



LARGE-ANGLE CORRELATIONS BETWEEN PROTONS 

Taking the calculated single-particle distributions f(f,p,t), we cannot make predictions for two- 
particle correlations. However, if these correlations are simply consequences of the conservation laws 
for momentum, energy, angular momentum, and particle number, then it is possible to obtain limited 
information on the two-particle correlation function, if one calculates in the ensemble method. It 
was shown (Bauer, 1987a) that two-proton correlation functions measured at large angles can be 
successfully reproduced by the BUU theory, provided that total momentum conservation is correctly 
taken into account. In another investigation Ardouin et al. (1988) found that the variation of the 
large angle correlation function with polar angle 9 is largely due to angular momentum effects. 
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CALCULATION OF CORRELATIONS AT SMALL RELATIVE MOMENTUM 



Even though it is impossible to produce two-particle correlation functions by only taking the cal- 
culated single particle distribution functions f(f,p,t), it is possible to approximately calculate the 
two-particle correlation function for small relative momenta, q = — P2), between the two emitted 
particles with momenta pf and p 2 . Here the interaction between the two particles has to be taken 
into account explicitly. 

To derive an expression for the two-particle correlation function, C(P, q), we assume that the final- 
state interaction between the two detected particles dominates, that final-state interactions with 
the emitting source and all remaining particles can be neglected, that the correlation functions are 
determined by the two-body density of states as corrected by the interactions between the two 
particles, and that the single particle phase space distribution function of emitted particles, g{p,x) 
varies slowly as a function of momentum p (i.e. g(p, x) ps g(p±q, x)). Then the theoretical expression 
for the two-particle correlation function can be written as (Koonin, 1977; Pratt, 1986; Gong et al, 
1991; Danielewicz and Schuck, 1992) 

C(P,q)=R(P,q) + l- U ^P^ 



n^oiLGry 

d 4 x 1 d i x 2 g{\P, x^g^P, x 2 ) 



q,n-f 2 +^l 



J d 4 Xl g{\P, Xl ) J d 4 x 2 g{\P, x 2 ) 



(7) 



where P = p\ + p 2 is the total momentum of the particle pair. x\ and x 2 are the space-time points 
of the emission of protons 1 and 2. IT is the single- and II12 is the two-particle emission probability. 

4>{q, r) is the relative wave function of the particle pair. The effect that gives rise to the HBT effect 
is the identical particle interference. In the absence of any other interaction the square of the two 
particle wave function is then simply given by 

|0(g,f)| 2 = 1 ±cos(2<ff) , (8) 

where the upper sign stands for bosons and the lower for fermions. 

In the presence of other interactions this result is modified. The relative wave function for Coulomb 
scattering is 

<l>c(q,r) = exp(-|7r?7) T(l + ir}) exp(iqz) 1 F 1 (-ir]\l\iq(r - z)) , (9) 
where 1F1 is the confluent hypergeometric series, 

. | | . ~ T(a + k)Y(b)z k . , 

and 77 = aZ\Z 2 m r j q. The two-pion (tt + tt + or 7r _ 7r _ ) relative wave function at small relative momen- 
tum is usually approximated as the symmetrized Coulomb scattering wave function 

|07™(g,r)| 2 = \\<pc{q,r) +<Pc{q,-r)\ 2 ■ (11) 



Due to symmetries, the wave function only depends on three independent variables which we choose 
to be q, r, and cos 6 = q- f/qr. The square of the two-pion relative wave function is displayed in fig. 
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Fig. 3: Absolute square of the two-pion (tt + tt + or n n ) relative wave 
function, {^^(q, r, cos6 l = 0.5)| 2 (From Bauer (1992b)). 



3 as a function of q and r for cos# = 0.5. For r —>■ 0, l^^l 2 is given by the the Gamov penetration 
factor 



|<Mg,0)| 2 = 2 ^ . (12) 

exp(27r?7) — 1 



For pion pairs of opposite charge (tt + tt ), the wave function does not have to be symmetrized, and 
it is given by the hydrogen-like Coulomb wave function with i](K + 7T~) = —ae 2 m r /q. 

For the two-proton relative wave function the strong interaction with the prominent 2 He-'resonance' 
cannot be neglected. To obtain the relative wave function in this case, we solve a radial Schrodinger 
equation with the Coulomb and the modified Reid soft core potential. The two-proton relative wave 
function is shown in fig. 4. One can observe the peak at q ~ 20 MeV/c due to the 2 He- 'resonance'. 
Due to the effect of antisymmetrization and due to the Coulomb interaction, \(j) pp \ — > as r — > 0. 

Equation |7| requires only the two-particle relative wave function and the smg/e-particle phase space 
distribution function. Under the assumptions stated above it is thus possible to generate two-particle 
correlation functions for small relative momenta from a theory which only predicts one-particle 
distribution functions. 

The resulting two-particle correlation function contains information on the space-time extension of 
the emitting source. To see this, we rewrite Equation as 

C(P,q) = J d 3 rF p (r) |0(g,r")| 2 . (13) 
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Fig. 4: Absolute square of the two-proton relative wave function, 
|0 pp (g,r,cos0 = O.5)| 2 . (From Bauer (1992b)) 



Here f = fi — f 2 is the relative coordinate of the two emitted particles, and the function Fp(f) is 
defined as 

/ d 3 R f{\P, R+lf, t>) f (±P, R- if, i>) 

*>(*0 = 1 — ~ 2 , (14) 

(J#r/(|P,f,*>)J 

where R = \{r\ + f 2 ) is the center-of-mass coordinate of the two particles, and the Wigner function 
f(p, f, t>) is the phase space distribution of particles with momentum p and position fat some time 
t> after both particles have been emitted: 

rt> 

f(P,r,t > )= dtg(p,r- p(t> -t)/m,t) . (15) 



For a given momentum P, the correlation function has three degrees of freedom, q, which are a 
function of Fp{f). Therefore correlation function measurements should allow the extraction of Fp(f), 

— * 

the normalized probability of two protons with the same momentum P/2 being separated by f. In 
this sense we are in principle able to extract information on the coordinate space extension of the 
emitting source at the emission time of the protons. 

One may also use correlation function measurements to test various theoretical models capable 
of predicting g(p,f,t) and thus making specific predictions about the correlation functions. This 
approach is more realistic in its goals, because a full six-dimensional determination of C(P, q) is very 
difficult in practice. 
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Fig. 5: Two-proton correlation functions calculated with a source param- 
eterization according to Equation [TJ| The solid lines represent 
the complete calculations including the effects of quantum statis- 
tics and of the Coulomb and strong interaction. The dashed lines 
represent the case of only Coulomb interaction, and the dotted 
line is for Coulomb interaction plus the effect of the Fermi-Dirac 
statistics for the two protons. (From Gong et al. (1991)) 



SENSITIVITY OF THE TWO-PROTON CORRELATION FUNCTION 

It is instructive to examine the sensitivity of the two-particle correlation function to different com- 
ponents of the two-particle interaction. We perform such a study for the two-proton correlation 
function. To do this, we use a simple zero-lifetime Gaussian source parameterization 

9o(p,r,t) = po exp(-r 2 /rl)5(t-t ) , (16) 

where r$ is the radius of the source. Figure 5 illustrates the effects of the different contributions 
to the two-proton final state interaction for sources of different radii. The Coulomb interaction 
dominates the shape of the correlation function for very large source radii. For r < 20 fm, the 
correlation function becomes increasingly sensitive to the effects of antisymmetrization and the strong 
interaction. The strong interaction has the dominant effect for source radii around ro = 2.5 fm due 
to the prominent 2 He resonance. 

It should be pointed out at this point, however, that in a realistic calculation the use of a simple 
zero life-time Gaussian source parametrization is not sufficient. Instead, calculations containing the 
full time and momentum dependence of the emitting source are needed. 
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Fig. 6: Two-proton correlation function for the reaction 14 N + 27 A1 at 
E/A = 75 MeV, as predicted by calculations based on the BUU 
theory (lines) and as experimentally measured (plot symbols). 
(From Gong et al. (1991)) 

CALCULATION OF TWO-PROTON CORRELATION FUNCTIONS 

We perform calculations of the single particle phase space distribution function / (r , p, t) by numer- 
ically solving Equation |[ These single particle distributions are then inserted into Equation [7] to 
generate the two-particle correlation function at small relative momentum. 

In Figure 6, we compare our calculations of the two-proton correlation function to experimental data 
for the system 14 N + 27 A1 at a beam energy of E/A = 75 MeV. Since the experimental data were 
not triggered on impact parameter, we have to also integrate our calculations over impact parameter 
with the proper weighting factors. (Details of this impact parameter averaging can be found in the 
Appendix of Gong et al. (1991).) The experimental and theoretical correlation functions are shown 
as a function of the relative momentum q for three different gates on the total pair momentum 
\P\ = |pi + P*2 1 - F° r comparison, the beam momentum per nucleon is Pbeam ~ 375 MeV/c. 

We show two different calculations. The solid line represents the full BUU calculations with the 
full nucleon-nucleon cross sections. The dotted line is the result of the BUU calculation with a 
reduced in-medium cross section, a = \a nni where a nn is the free space elementary nucleon-nucleon 
cross section. In both cases the medium corrections due to the Pauli-principle for the final nucleon 
scattering states are, of course, taken into account. 

It is clear from this figure that sizeable differences between the two calculations with different as- 
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Fig. 7: Sensitivity of the two-proton correlation function to the in- 
medium nucleon-nucleon cross section, a, and to the nuclear 
compressibility. (After Gong et al. (1990)). 

sumptions on the in-medium nucleon-nucleon cross section exist. In fact, variations of the cross 
section by only 10% result in differences between the calculated correlation functions which should 
be experimentally measurable. From this one has to conclude that two-particle correlation functions 
are very sensitive probes for the collisional dynamics of intermediate energy heavy ion collisions. 

In fig. 7 we display the sensitivity of the two-proton correlation function to the in-medium nucleon- 
nucleon cross section, a, and to the value of the nuclear matter compressibility, k. We show the 
results of our calculations for a stiff nuclear equation of state (k = 380 MeV) and different values of 
the in-medium nucleon-nucleon cross section. 

From our theoretical results and comparisons to experimental data we conclude that the value of the 
in-medium cross section (aside from correction due to the final state 'Pauli blocking') is very close to 
the experimentally measured (energy dependent) free value. This is in agreement with information 
we extracted from our investigation of the disappearance of nuclear collective flow (Krofcheck et al, 
1989, 1992; Ogilvie et al, 1990; Bauer 1992a). 

We also varied the compressibility of nuclear matter, which enters the calculations through the density 
dependence of the mean field potential U . In fig. 7, we also show a calculation for a soft equation of 
state (k = 200 MeV). Here we find, however, only a weak sensitivity of our results on this parameter. 
This is expected, because for the relatively light system and low beam energy considered here only 
moderate maximum values of the nuclear density are achieved. 
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Fig. 8: Comparison of theoretical calculations (histogram, Li and Bauer, 
1991a) and experimental data (circles Odyniec et al., 1988) for 
single pion kinetic energy spectra from the reaction La+La at 
E/A = 1.35 GeV. 



CALCULATION OF TWO-PION CORRELATION FUNCTIONS 



Pluta et al. (1992) have observed that when one computes the ratio 

a 2 (TT + ,TT-) 

<7i(7T+)<7i(7T-) 1 ' 

as a function of the invariant mass of the pion pair for 1.6 GeV protons incident on heavy target 
nuclei, there is a striking enhancement close to M inv (7r + , n~) = 2M W . 

We have investigated the possibility that the enhancement of the correlation function close to an 
invariant mass of 2m w could be caused by the effect of focussing of the two correlated pions into 
the same hemisphere due to the presence of the target spectator matter, which causes some of the 
produced pions to be rescattered and absorbed. To do this, we performed transport calculations 
with an extended BUU code, which had previously been successfully used to calculate pion spectra 
in heavy ion collisions (Li et al. (1991)), and which has been used to show that the observed pion 
collective flow (Gosset et al. (1991)) was caused by nuclear shadowing effects (Li et al. (1991b)). 

In this extended BUU code we propagate protons, neutrons, Delta and N* resonances, and pions. 
For their interaction cross sections, we use the experimentally measured free hadron-hadron cross 
sections where available. In the cases where the elementary cross sections are not measured we 
employ isospin symmetry and/or detailed balance to obtain the unknown quantities. The process of 
greatest interest here is the creation and reabsorption of pions, which is dominated by the two-step 
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Fig. 9: Two pion (ir + , ir~) correlation function as a function of the invari- 
ant mass of the pion pair for p + C, Pb reaction at an incident 
energy of 1.6 GeV. The data are from Pluta et al. (1992) and 
are represented by the squares. Our calculations (Klakow et al. 
(1992)) are represented by the histograms. 



Figure 8 (Li and Bauer, 1991a) shows the single pion spectra calculated in this model. The reaction 
is La+La — > 7r~+X at E/A is 1.35 GeV. The data of Odyniec et al. (1988) are represented by the 
circles, and the results of the calculation are shown by the histogram. Experiment and theory were 
both gated for central collisions (b < 2.8 fm), and the detection angle for the pions was 90° ± 30° in 
the center of mass of target and projectile. More recent calculations (Li et al, in preparation) use 
improved parametrizations for the elementary cross sections and shown an even better agreement 
with experimental data. 

To generate correlation functions which are comparable to experiment, we calculate the correlated 
pairs from events with fixed reaction plane 



process 



N + N^N + A^N + N + tt. 



(18) 




(19) 



whereas we compute the uncorrelated background events from 




(20) 
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where E { and pi are the energy and momentum of pion i in the some frame of reference (we chose 
the lab frame). The operator IZ^ rotates a vector p by a (random) angle about the beam axis (the 
z- axis in our case): 

cos — sin \ ( p x 
T^4>P — | sin0 cos0 \ Py \ (21) 
1 / \ Pz 



The correlation function is then computed from 

iV((M inv ^) 



= TTTT^M , (22) 



With this technique we approximate the experimental fact that for the generation of the background 
events one is forced to average over random relative orientations of the reaction planes for the two 
uncorrelated pions in every pair, whereas the correlated pairs are all sampled with fixed reaction plane 
of the pair and therefore contain the effects of the target matter shadowing at finite impact parameter. 
For a more detailed description of the numerical procedure and the effect of pion absorption on two 
pion correlation functions see Klakow et al. (1992). 

In fig. 9, we compare the results of our calculations with the experimental data for the C and Pb 
targets. We find that we can reproduce the Pb target data by the calculations, whereas we slightly 
overpredict the correlation function at small M inv for the C target. 

The main result of this investigation is thus that the effect of shadowing due to pion absorption has 
a large effect on the extracted two-pion correlation function, particularly in asymmetric projectile- 
target systems. As shown this kinematic effect leads to an enhancement of the two-pion correlation 
function at small M inv and can lead to wrong extracted source sizes if not included. 



CONCLUSIONS 

The calculation of two-particle correlation functions at small relative momentum on the basis of 
one-body transport theories is feasible by using the convolution techniques described above. Thus 
intensity interferometry is a powerful tool to test nuclear transport theories and to investigate nuclear 
dynamics. 

Comparisons with experiment show that the BUU transport theory is able to reproduce detailed 
features of the experimentally measured two-proton correlation functions. These features include the 
pair momentum dependence of the peak due to the 2 He-'resonance', the effect of source deformation, 
and the lifetime effect on the correlation function. 

We have shown that the theoretically obtained two-proton correlation functions are sensitive to the 
value of the in-medium nucleon-nucleon cross section. At higher beam energies and for large systems 
we also expect a sensitivity of the results on the compressibility of nuclear matter. Thus nuclear 
intensity interferometry is a useful tool to investigate the nuclear transport properties, and it should 
also enable us to conduct further studies of the nuclear equation of state. 

From a theoretical standpoint it is clearly desirable to compare to impact parameter resolved exper- 
imental data, which will further increase the sensitivity of two-particle correlation functions at small 
relative momentum to the effects discussed above. These studies are currently in progress. 
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For the two-pion correlation functions we have shown that it is essential to included the effects of 
pion rescattering and absorption properly into the calculations. A source characterization without 
taking these effects into account will lead to wrong results for quantities like the size of the emitting 
source. 

For the next few years, we anticipate that the focus of HBT type of studies in heavy ion collisions 
will also shift to the energy range of E/A pa 1 GeV. At this energy range we expect much higher 
compression (p/po ~ 2 — 4) of the colliding nuclei than at the beam energy range considered here 
(E/A pa 100 MeV, and p/po ~ 1 — 3 — 1.5). Then two-particle correlations should be able to give 
us additional information on the nuclear compressibility and thus provide valuable insight into the 
nuclear equation of state. Work in this direction is alrady in progress (Mader and Bauer, 1992). 
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